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Abstract 

Chiral effective field theories have been used with success in the study of nuclear 
structure. It is of interest to systematically improve these energy functionals (par- 
ticularly that of quantum hadrodynamics) through the inclusion of many-body cor- 
relations. One possible source of improvement is the loop expansion. Using the 
techniques of Infrared Regularization, the short-range, local dynamics at each or- 
der in the loops is absorbed into the parameterization of the underlying effective 
lagrangian. The remaining nonlocal, exchange correlations must be calculated ex- 
plicitly. Given that the interactions of quantum hadrodynamics are relatively soft, 
the loop expansion may be manageable or even perturbative in nuclear matter. 
This work investigates the role played by the three-loop contributions to the loop 
expansion for quantum hadrodynamics. 
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1 Introduction 



Density Functional Theory (DFT) is a powerful technique originally developed 
for use in condensed matter physics pl2l[3] that has been successfully adapted 
to nuclear physics [4f5]l6f7f8|9] . DFT states that the ground-state expectation 
value of any observable is a unique functional of the exact ground-state den- 
sity; moreover, if the expectation value of the hamiltonian is considered as a 
functional of the density, then the exact ground-state density can be deter- 
mined by minimizing the energy functional p!|3] . Furthermore, DFT allows 
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one to replace the quantum many-body equations by a series of single-particle 
equations with local, classical fields that reproduce certain observables exactly 
(energy, scalar and vector densities, and chemical potential) |2li9iil0j . Thus, the 
problem is reduced to determining the exact ground-state energy functional. 
However, this is impossible in practice. As a result, a number of approximate 
energy functionals have been developed; one such theory is based on quantum 
hadrodynamics (QHD). 

QHD is a low-energy theory of the strong interaction [7f5f51[TT|12lll3lll4|15j . 
Here the hadron, and not the quark, is the observed degree of freedom (due to 
confinement at this energy scale). QHD models the nuclear force as a exchange 
of mesons between nucleons. Isoscalar scalar (cr) and vector {u) mesons repre- 
sent a medium-range attraction and a short-range repulsion respectively. The 
pion is also included to take chiral symmetry into account. DFT allows one to 
replace the quantum meson fields with their classical equivalents (Kohn-Sham 
potentials called mean fields); these mean fields, while large, are small com- 
pared to the chiral symmetry breaking scale. As a result, one can use these 
ratios as small parameters with which to expand the energy functional in a 
controlled fashion. Each term in this lagrangian is characterized by an unde- 
termined coefficient which is assumed to be natural, or order unity [TBPT] . The 
resulting lagrangian [7f8] has provided a method for predicting the properties 
of nuclei fT8][T9]l20ll2T]|22]|23l . 



However, the mean field theory of QHD is only an approximation to the exact 
energy functional. The question remains how this theory behaves when other 
many-body corrections are included to improve the energy functional (such 
as loops, rings, clustering, etc.). This work is part of an investigation into 
the effects of loops on this particular energy functional and is a continuation 
of [21I2S]- In these works, the many-body loop expansion was carried out to 
the two-loop order where, using the techniques of Infrared Regularization, all 
of the short-range, local dynamics were absorbed into the parameterization of 
the underlying effective lagrangian. What remained was the nonlocal exchange 
correlations, which were then explicitly calculated. In addition, the effect of 
this expansion on naturalness was investigated and the size of the two-loop 
exchange integrals was determined. Since the interactions of QHD are rela- 
tively soft, one might expect that, once Infrared Regularization is taken into 
account, the loop expansion may be asymptotic. It is the purpose of this work 
to investigate the effect of the three-loop contributions on this theory. 

The loop expansion for QHD is constructed in the usual manner [2^27ll28|29ll3Ulli 
The effective action is expanded around its classical value by grouping terms 
according to the number of quantum loops (or powers of h) in their corre- 
sponding diagrams. For the purposes of this work, we are interested only in 
the terms in the expansion at the three- loop level. All of the integrals that 
represent tadpole and disconnected diagrams cancel out in the effective action 
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and we are left with only the fully connected diagrams. In addition, those 
diagrams which are anomalous are discarded. For the cases considered here, 
there are eleven integrals of interest. Ten of these integrals have four factors 
of the nucleon propagator and two meson propagators (either scalar, vector, 
or pion) each. There is one additional three-loop diagram with two factors of 
both the baryon and the pion propagator. In this work, no nonlinearities in the 
isoscalar mesons were included in the effective lagrangian, as in [23]. In [25], 
the effect of these nonlinearities was explored at the two-loop level. As we are 
interested only in the general effects of the loop expansion at third order (and 
not an improved equation of state), these nonlinearities are not retained. The 
effect of the inclusion of nonlinear meson self-interactions on the three-loop 
integrals is left for future work. 

Why a loop expansion? The loop expansion is a simple and well-developed 
expansion scheme in powers of h that is derived from the path integral. The 
mean meson fields are included non-perturbatively and the correlations are 
included perturbatively. Therefore, one can analyze the many-body effects or- 
der by order. The loop expansion has the advantage that it is fairly easy 
to separate the local and nonlocal dynamics and analyze their structures 
|24f25|33|34|35|36ll37] . Consideration of the three-loop integrals will help us 
determine if this is a good expansion scheme or not for QHD. 

The diagrams at the three-loop level can be written as either ladders or crossed 
ladders (with the exception of the one additional pion diagram). If one consid- 
ers the Bethe-Saltpeter equation in the ladder approximation, one can replace 
the two nucleon propagators in the ladder with an approximate two-body 
propagator |29|38|39|40|41II42] . This two-body propagator contains the nonlo- 
cal physics from the two one-body propagators and provides the proper cuts 
to maintain analyticity. It describes the propagation of two nucleons above the 
Fermi surface. As there is some ambiguity in the definition of this propagator, 
we will consider three different versions: Blankenbecler-Sugar, Thompson, and 
Erkelenz-Hohnde [29lliO]liTE2] . 

The remaining pion (football) diagram arises from the inclusion of a two- 
pion vertex. Isospin and parity considerations restrict the two pion vertex to 
diagrams with the pions and nucleons only, at least at the three-loop level. 
Because this vertex is antisymmetric, only a diagram that has two of these 
vertices will survive. Fortunately, this diagram is analogous to the ladders and 
crossed ladders; here we use a two-body pion propagator (essentially the same 
as the baryon case but without the projection operators and theta functions) 
. Note that this two-pion exchange is not equivalent to scalar [a) exchange 
(which is an effective field representing any interaction in the isoscalar-scalar 
channel). Any potential redundancy is eliminated by the parameterization. 

The remaining nucleon propagators can be separated into two components: 
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Feynman and Density parts [26f29] . This is accomplished by taking into ac- 
count the proper pole structure of the propagator. The Feynman part describes 
the propagation of a baryon or an antibaryon; the Density portion describes 
on-shell propagation of a nucleon while taking the exclusion principle into ac- 
count. Thus, one can separate the local and nonlocal dynamics, as in p^l[25] 
(see also Infrared Regularization discussions in j33|34|35ll36ll37j ). All of the 
vacuum loops are parameterized; only the boson exchange between valence 
nucleons is calculated. 

It is also of interest to investigate how these loop integrals fit into the power 
counting scheme that was developed for the underlying QHD lagrangian 
[71I51I14II15] . It was shown in previous work that the two-loop integrals are 
about third order in the power counting [2H[25] . In this work, we will deter- 
mine the size of the three-loop integrals and how they relate to both the one- 
and two-loop levels. The tricky business of finding the appropriate expansion 
parameter is left for future work. 

The loop expansion has been studied in the context of other effective field the- 
ories, particularly Chiral Perturbation Theory (ChPT) |33f34ll35f36ll37ll43f44j . 
The techniques of Infrared Regularization were developed to separate out and 
absorb the local portions of the loop expansion in ChPT [53|34ll35ll36f37j . 
More recently, investigations have been conducted to determine the effect of 
multi-loop correlations in ChPT and to discover the appropriate expansion 
parameter [13|1^ . 

In this paper, the many-body loop expansion for QHD is investigated at the 
three-loop level. Infrared Regularization allows us to parameterize the local 
dynamics. The nonlocal exchange correlations are solved in the ladder approx- 
imation to the Bethe-Saltpeter equation for the lowest-order truncation in the 
underlying lagrangian. The effects of a high momentum cutoff are explored. 
The naturalness of the lowest level couplings is confirmed at the three-loop 
level. In addition, the size of the integrals is discussed and how they fit into 
the power counting scheme is investigated. 
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In this section, we present the corrections to QHD arising from the three-loop 
contributions to the loop expansion. The background for the loop expansion 
in QHD is presented in previous works pi|25ll26] . The lagrangian used here 
follows from the chirally invariant lagrangian in [7.8], or 



C=-%Ij [7^ [df, - igvVf,) - igAjf^JbCifi + hf^Vf, + (M - gs(p)] ip 
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where vr = ■ Ta. Here ip are the fermion fields and (p, V^, and Ha are 
the meson fields (isoscalar-scalar, isoscalar- vector, and isovector-pseudoscalar, 
respectively). The heavy meson fields are also chiral scalars. Note that in this 
work, the conventions of [9] are used. As in pH[25] , we are not making a chiral 
expansion in powers of the pion mass and include it for kinematical purposes 
only. The pion-nucleon interactions retained above can be expanded in the 
following fashion (with ^ = exp{— iTr/ZTr}) [7ll8] : 



^, = '1 + id,(!) = k, d,iL] + -- - , (2) 

«M = - ^ - id,e) = -jd,n + ... . (3) 

It is necessary to retain only the first term in and at the three-loop 
level. Note that the first term in f ^ is a two-pion-nucleon vertex that does not 
contribute at the two-loop level. 

The generating functional is defined in the usual way by 



Z[j,J^]=exp{tW[j,J^]/h} 

= Ar-i| D{^)D{ij)D{<j))D{V^)D{rra) 

xexpl^^J d'^x[C{x) + j{x)(t){x) + J^{x)V^{x)]^ , (4) 

where Af is the normalization factor (in effect, the vacuum subtraction), j{x) 
and J^{x) are the external sources corresponding to the meson fields and 
V^, respectively, and the connected generating functional is W[j, J^]. 



2.1 Ladders and Crossed Ladders 



Consider the portion of the connected generating functional that describes the 
scalar-scalar three-loop contributions, which is 



3-55 '■ 



9s J J d'^xd^yd'^zd'^a 



'24 
6u{x) 





-i5 




-i6 






5u{y) 




6u{z) 





-i6 
6u{a) 
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X exp |-i J J (fxid!^X2i{xi)G h{xi - a;2)^(x2)| 



xexp|- 
-VEV . 



d'^xid^X2u{xi)^%{xi — X2)u{x2] 



sources=Q 



(5) 



We then perform the variational derivatives and drop the terms with the 
tadpole diagrams (these diagrams will cancel when we consider the effective 
action). Three of the remaining terms are disconnected diagrams (the ones 
that are essentially the scalar two- loop integral squared); these are discarded 
for the same reason. Next, we perform Fourier transforms over all the terms. 
The three delta functions that arise let us eliminate three of the six momentum 
integrals. Rearranging and combining terms, we get (after transforming to the 
energy density and suppressing the h) 



~^Ti J (27r)4 (277)4 (27r)4 

X {^s{k)^s{k) tr [Gniq - k)GH{q)] tr [GH{q')GH{q + k)] 
-As{k)As{q' - k) tr [Gniq' + q - k)GH{q)GH{q + k)GH{q + q)] 
-2As{k)As{q - q') tr [Gniq - k)GH{q)GH{q')GH{q)]} ■ 

(6) 



Note that the propagators used here are defined in [21]. These integrals cor- 
respond to the three diagrams in Fig. [H The term that produces the third 
diagram in Fig. [1] drops out at zero temperature because it produces a particle- 
hole pair of equal momentum at the Fermi surface, which violates the Pauli 
exclusion principle (one can see this more clearly if the meson propagators are 
shrunk to contact interactions as in Fig. [21) [45] . The first and second diagrams 
in Fig. [1] are in fact ladder and crossed ladder diagrams respectively, as will 
be seen later. The process by which we have arrived at the above integrals 
can be repeated for the various possible combinations of meson propagators. 
This will yield a ladder and a crossed ladder for each of following combina- 
tions: scalar-scalar, scalar-vector, vector-vector, scalar-pion, vector-pion, and 
pion-pion. However, the scalar-pion and vector-pion ladders do not conserve 
parity and therefore vanish. This leaves a total of ten integrals (4 ladders and 
6 crossed ladders) which are shown in Fig. [31 Note that the above analysis 
has taken into consideration only the one-pion vertex [the term proportional 
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to in Eq. (^]. 




Fig. 1. Scalar-scalar three- loop diagrams defined in Eq. Q. Here the double line 
represents the baryon propagator and the dashed line represents the scalar meson 
propagator. 



Fig. 2. Scalar-scalar three-loop diagrams in which the meson interactions were 
shrunken to a point. The first and second diagrams in Fig. [T] appear the same 
here although there are two ways to do the trace(s). 

Now, we can write the connected generating functional which contains all the 
ladders (L) and crossed ladders (XL) at the three-loop level. Proceeding as 
before, the remaining terms are (where dimensional regularization was used 
to make the substitution Gh G* [2lll25ll26]l29l| ) 



S,iL + XL) ' ' 



(27r)4 (27r)4 (27r)4 
X i^-z^As{k)As{k) tr [G*{q - k)G*{q)] 

xtr [G*{q')G*{q' + k)] 

+t^Asik)Asiq' - k) tr [G*iq' + q - k)G*{q) 

X G*{q + k)G*{q' + q)] 

+z^-^Asik)V%ik) tr [G*{q - kh,G*{q)] 

xtr [G*iq'h,G*iq' + k)] 
_,MA5(fc)I)J^(g' - k) tr [G*iq' + q - k)^,G\q) 

X G*(g + A;)7,G*(g' + g)] 
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-^^V%{k)Vl{k) tr [7,G*(g - k)j^G*{q)] 
xtr [7,G*(g')7.G*(g' + A;)] 

+z^P;,(fc)Pl(g' - A;) tr [^,G*{q' + q - k)^.G*{q) 



X7,G*(g + A;)7^G*(g' + g)] 



,2 „2 



+i||fAc,(fc)A;^(g'-A:) 

J TT 



X tr 



G*(g' + g-A:)(/- }i)^i,G\q) 



xG*iq + k)i^'- ^)|75G*(g' + g) 



X tr 



If^Gniq' + q 



-l5GH{q) 



Xl.GH{q+k){^'- f^)p,GH{q' + q) 



—I 



A^^ik)Al\k)tT 



r^^G*iq-k) ^75yG*(g) 



X tr 



9\ 



)^^G\q') fi^,jG*iq' + k) 



+z^A;^(fc)At'(g'-A:) 



X tr 



^75^G*(g' + g-fc)(/- fiHjG*{q) 



X l^l,^G*{q + k)U'- ]^)^J-iG\q' + q) 



(7) 



In Fig. [3], the diagrams corresponding to the integrals in Eq. (J?!) are shown, 
in the same order. 



To calculate these ladders and crossed ladders, we must first consider the 
Bethe-Saltpeter equation. Two-body correlations can be introduced through 
an effective interaction F, which is the solution to the Bethe-Saltpeter equation 
in the nuclear medium [38|39ll47j 



T = K + J KGGT 



(8) 



where K is full two-body scattering kernel and G is the fully interacting baryon 
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Fig. 3. Surviving three-loop diagrams redrawn to show exphcitly that they are in fact 
ladder and crossed ladder diagrams (excluding diagrams with the two-pion vertex). 
Here the double lines represent the baryon and the dashed, wavy, and dotted lines 
correspond to the scalar (o"), vector (lj), and pion respectively. Notice that there 
are no scalar-pion and vector-pion ladders because they do not conserve parity. 

propagator. In practice, the full kernel cannot be written in closed form; we 
replace it with the usual ladder approximation, K ^ V. Here V is just a 
one-meson-exchange potential. Then Eq. ([8]) is rewritten as two equations, or 



where g is an approximate two-body propagator and U is a quasipotential. 




(10) 



(9) 
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For the three-loop case, these equations become 



T = U + I UgV , 

U = V + jv {GG -g)U . 



(11) 



(12) 



If the two-body propagator [g) is a good approximation for the product of 
the two interacting baryon propagators, then Eq. f|T2l) imphes U and Eq. 
ffTTl) becomes 



Thus we can replace two of the baryon propagators in each of the ladder inte- 
grals in Eq. ([7]) with a two-body propagator. Specifically, we want to replace 
the two propagators that form the interior loop in the ladder diagrams in Fig. 
[3l this turns out to be the loop over the momentum k. For the case of the 
crossed ladders, we will also insert the two-body propagator into the k loop; 
here however, we will also need to perform a Feirz transformation to properly 
antisymmetrize the integral [15] . 

We replace the two baryon propagators in the k loop with the Blankenbecler- 
Sugar two-particle propagator (assuming G*G* ~ g) [29ll38f 3911^ . 



where the total momentum coming into the k loop is P^, P* = + 2S^, and 
s* = P*^ and is assumed to be constant with respect to the loop momentum 
[39|l^ . The is included to eliminate the self energy contributions to the 
momentum. As in the two-loop case, we have 




(13) 



gik,P*\B) = -n6ik,) 



Q{k,P*\B)K^l\\P* + k)K^l\\P* - k) 



(14) 



E*{k) E*^{k) - Is* -ie 




(15) 



We also define the following projection operators 




(16) 



A^^\^P* -k)=tj. (Ip* -k)- i^E*i}-P* -k)-M* . 



(17) 
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In the nuclear matter limit, we write the Pauli exclusion operators as 



Q{KP*\B) 



l-d[kp- \-P + k 



kp- \-P-k\ 



Here P* = P (this arises from a relativistic generalization of the "reference 
spectrum" approximation that is used in nonrelativistic Brueckner calcula- 
tions Eq. (ITH]) ensures that the two interacting baryon propagators 
in the k loop are above the Fermi surface. The propagator in Eq. fll4p is an 
approximation for the two baryon propagators that provides the proper cuts 
to ensure analyticity and removes the local effects that Infrared Regulariza- 
tion eliminates. In addition, other versions of the two-body propagator can 
be used; the alternate versions considered here are the Thompson two-body 
propagator [25|1T] 



g{k,P*\B) 



-nd[k4) 

2 E*{k)^/^ 



(1), 



IP* + k)A^+\lP* - k) 



E*{k) + -ie 



(19) 



and the Erkelenz-Holinde two-body propagator [29IH2] 



g{KP*\B) 



1 



-T,b{k4^-^*-E*{k)) 



X 



Q(k, P*\B)t^l\\P* + k)t^l\\P* - k) 



E*{k) E*^{k) - is* - ie 



(20) 



Now we are ready to consider the scalar-scalar ladder diagram (the first dia- 
gram in Fig. [3]); the box in the center is the loop over k. The ladder diagram 
can be written using the two-body propagator in Eq. (fT4l) (with P^ = q'n + g^) 



3-SS 



9s^ 



X 



d^k d^q d^q' 
(27r)4 (27r)4 (27r) 



-\s{k)/\s{k)Q{k,P*\B) 



5{ki) tr 


'G*{q)kf{\P*-k) 


tr 


G*{q')A^^\lP* + k) 


E*{k) 


E*'^{k) - \s* - ie 





(21) 



Now we let G* G*p + [261129] . The terms with one or more factor of G*p 
will not have enough delta functions to eliminate all the frequency integrals 
and are therefore expressible as a polynomial of terms that are already present 
in the underlying lagrangian. These terms are just absorbed and do not need 
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to be calculated. That leaves us with only the nonlocal portion. We can align 
the k momentum along the z-axis; this allows us to integrate the angular 
portions of the k integral. Now, working out the traces (which are over both 
spin and isospin) and writing out the angular integrals, we get 



kp 



^3-ss{L) 



1287r8 J EHq) J E*(q') 

^ -1 



gl 7 7 \qfd\q'^ 



d{cos9gk) / d{cos9g'k 



-1 



2tt 2n 

X J d(j)gk J d(f)g'k 





\k\'^d\k\ 
E*{k) 



Q{k,P*\B) 



X 



X 



E*{q)E*{]^P* -k) + M*^ - q- Qp* - k 
E*{q')E*{^P* + k) + M*^ - q' ■ Qp* + k 



/ \k\' + ml 



\k\^ + M*"^ - -s* - ie 



(22) 



The scalar-scalar crossed ladder diagram (the second diagram in Fig. [3]) with 
the two-body propagator in Eq. (fT4l) is (with = q'fj^ + 2q^) 



d^k d^q d^q' 
(27r)4 (27r)4 (27r) 



-As{k)As{q'-k)Q{k,P*\B) 



X- 



5(A;4)tr 


Aj^(iP*- 


k)G*{q)A^^\lP* + k)G*{q' + q) 


E*{k) 


E*'^{k) - \s* -ie 





(23) 



Now we let G* G*, 



G*j~i [26|29] . As before, the terms with one or more fac- 



tor of G*p will not have enough delta functions to eliminate all the frequency 
integrals and therefore are expressible as a polynomial of terms that are al- 
ready present in the underlying lagrangian. These terms are just absorbed and 
do not need to be calculated. That leaves us with only the nonlocal portion. 
We can align the k momentum along the z-axis; this allows us to integrate the 
angular portions of the k integral. However, the trace is not as straightforward 
as before. In order to use the two-body propagator, we must also invoke the 
Fierz identity. This essentially antisymmetrizes the ladder diagram. The Fierz 
identity is [i6] 



^ X P: 



AP X / + A-i^B X 75 + A-i^B X 7„ 
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(24) 



where A and B represent either 7^ or / and the x represents any matrices 
that come between A and B. Substituting this into the scalar-scalar crossed 
ladder integral, we get 



9i fmoi 



2%8 7 E*(q) 

-1 



d{cos6qk) / d{cos6, 



2n 2n 



'Jq'k 



X 



X 



X 



\q'?d\q'\ 

E*{q' + q) 



0{k,-\q' + q^)j\^Qik,nB) 



ip* - Pj ■ [q' + g) - E*(ip* - k)E\q' + q) - M 
ip* + fcj • g - E*{^P* + k)E*{q) - M 



* 2 



*2 



g'l-l^l ~[E*{q' + q)-E*{q)Y + m's 



X 



-s — te 
4 



(25) 



The other ladders and crossed ladders are handled in a similar fashion. The 
only difference is the traces are more complicated and the coupling constants 
change. In addition, one can substitute Eqs. (fT9l) and (l20l) for Eq. (JT^ . 



2.2 Two-pion Vertex 



In this subsection, we consider the contributions to the loop expansion at 
the three-loop level involving the two-pion vertex. These terms arise from the 
second pion-nucleon coupling in Eq. ([T]). There are two possible combinations 
that conserve isopsin and parity (and are not redundant); this corresponds to 
the following portion of the connected generating functional: 



3 — TTTT ■ 



2 2/4 



d'^xd'^yd'^zd'^z' 
-i5 



dUx) 



n 
2 



-i5 



Kk{z) 



gz 



aa' 



46 



2 ' 2 



6zz' 



11 
2 



/3/3' 



77 
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X 



iS 





-iS 




iS 




-iS 




iS 




- a 




a' 


my)\ 




my)\ 


13' 




7 - 



-iS 



.h' 1 



2 4/^ 
X 17m 



J J d^xd^x'd^ydy 
-id 



-i6 

-id 



-i6 



X 



X exp 



i6 





-i6 




i6 




- a 




a' 


my)\ 


p - 



.2 ' 2 

Tk n 

2 ' 2. 
-i6 



6yyl I 



d'^Xid'^X2^{xi)GHixi — X2)^{X2 

exp 1^ y" J d'^xid'^X2Cc{xi)A'^{xi - X2)Cd{x2) 



X 

-VEV 



sources=0 



(26) 



These terms are represented by the diagrams in Fig. |H After working out the 
variational derivatives, taking the Fourier transforms, and using all the delta 
functions, the term corresponding to the second diagram in Fig. H] drops out. 
As the two-pion vertex is antisymmetric, only a combination that includes 
two of these vertices can survive. The remaining (football) diagram is repre- 
sented by the following integral (after transforming to the energy density and 
suppressing the h) 

/r^\ ■ ^ f f d'^q d^q' . ,,, . , 

^3-..(i^) = y y (2^ (2^ (s^^'^^^)^-^'^ ) 

xTr{(/- ^)G*(g) ^G*(g-A;-g')} , (27) 

where Tr is summed only over the spin (the isospin sum has already been 
done). Again, dimensional regularization was used to make the substitution 




Fig. 4. Three-loops diagrams involving the two-pion vertex. 
We now proceed in a manner similar to the method used in section 12. 1[ Here 
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we need to introduce a two-pion propagator (A^rA^ ~ A^rTr), 



AUk,P*) = -7l6(h) T 

^ '^Em Efik) 



1 „* 

4'* 



m 



where 

E:{k) 



fcp + ml 



(29) 



The two-pion propagator follows from the two-body Blankenbecler-Sugar nu- 
cleon propagator without the projection operators and the Heavyside step 
functions pHIHU] . Thompson and Erkelenz-Holinde forms can also be used. 
This allows us to rewrite the integral as (with P^ = q'^) 



^3-Tnr{E) 



d^k d'^q dS' 
16/4 7 J (2^(2^(M 



-Mk, 



xTr[(/-2^)G*(g)^G*(g' + g)] 



/\E:{k) 



Ef{k)--s*~ie 



(30) 



G}) [26|29] . As before, the terms with one or more fac- 



Now we let G* 

tor of G*p will not have enough delta functions to eliminate all the frequency 
integrals and therefore are expressible as a polynomial of terms that are al- 
ready present in the underlying lagrangian. These terms are just absorbed and 
do not need to be calculated. That leaves us with only the nonlocal portion. 
We can align the k momentum along the z-axis; this allows us to integrate the 
angular portions of the k integral. Now, working out the trace (which is only 
over spin) and writing out the angular integrals, we get 



^ ^ ^ d{cOs9qk) j d{cOs9qik) j d(j)gk J t 
-1 

(lfd\^'\ nfu , ^\ f \k\^d\k\ 



X [P {2E*{q) [E*{q' + q) - E*{q)] - 2q' ■ q} 
+q ■ k <Aq ■ k — 2q' ■ q + 4q' ■ k — q' 



+ [E*{q' + q)- E*{q)] [E*{q' + q) + E*iq)] }) 
/{Ef{k)-y + \[E*{q' + q)-E*{q)f 



(31) 
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Fig. 5. The total three-loop energy density as a function of cutoff. Three versions of 
the two-body propagator were investigated: Blankenbecler-Sugar (BS), Thompson 
(T), and Erkelenz-Holinde (EH). Here the MOA parameter set was used |24] . 

3 Discussion 

In this section, we discuss the results of the three-loop calculations for QHD. 
The level of truncation retained is sufficient to study the general effects of 
the three- loop energy (as can be seen at the two- loop order in No 
nonlinear meson self-interactions were included in the underlying lagrangian. 
We leave the inclusion of these nonlinear effects for future work. In addition, it 
was shown implicitly that the local, short-range dynamics were absorbed into 
the parameterization of the lagrangian (as it is expressible as a series of terms 
that in principle already exist). What remains are the nonlocal, long-range 
correlations, which are calculated explicitly. 

Now we consider the numerical analysis of the three-loop energy. We use the 
parameter sets listed in Tabled! For the pion terms, we use g\ = 1.5876 and 
U = 93 MeV 0. In addition, both the Wl (one- loop level) [7j and MOA (two- 
loop level) [21] sets lead to equilibrium at kp = 1.3 fm~^. The mean meson 
fields are determined by extremizing the meson field equations and are used 
as input to the exchange integrals. Then the full three-loop energy density 
is extremized with respect to the meson fields. The two-loop integrals were 
calculated using Gaussian quadrature. 

A total of eleven integrals survive at the three-loop level: four ladders, six 
crossed ladders, and one football diagram. The three-loop integrals are also 
solved using Gaussian quadrature. However, the interior loop (the one over 
momentum k) in all eleven integrals has an upper limit of infinity. As a result, 
a cutoff must be introduced. This is expected, since the effective field theory 
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is applicable only up to some breakdown scale determined by non-Goldstone 
boson physics. The sum of the three-loop contributions to the energy density 
as a function of cutoff is shown in Fig. O The results in Fig. O are for three 
different versions of the two-body propagator [see Eqs. (HM . (HM . and (I2U]) ] 
and utilize the MOA parameter set [21]. The choice of the cutoff is critical. If 
the cutoff is too high, then some local physics gets included and the energy 
becomes large. If the cutoff is too small, some of the nonlocal dynamics is 
lost. A reasonable cutoff (600 < A < 700 MeV [TSElj ) yields a total three- 
loop energy of less than 17 MeV. For comparison, the total two-loop energy 
for the MOA set is 29.60 MeV. One can see that for a reasonable choice of 
cutoff, the energy decreases from the two- to three-loop order. 

It is clear from Fig. [5] that there is some uncertainty arising from the two- 
body propagators. To get an improved idea of how good an approximation 
each propagator is, Eqs. (ITT]) and (IT^ need to be solved self-consistently. This 
is left for future work. 

In Table |2l the individual contributions are shown using the MOA set for four 
choices of the cutoff (A = 500, 600, 700, and 800 MeV) and the Blankenbecler- 
Sugar propagator. The terms are quite large; however, there is always some 
cancellation between ladders and their corresponding crossed ladders. In ad- 
dition, there is large cancellation when the scalar-scalar, scalar-vector, and 
vector-vector terms are summed. 

The set MOA listed in Table [T] was fit at the two-loop level to nuclear equi- 
librium {£ / Pb — M = —16.1 MeV and kp = 1.30 fm"^) by adjusting gs and 
gv using a downhill simplex method to minimize a least-squares fit (with re- 
spective weights of 0.0015 and 0.002) [24j. The sets MOBl and M0B2 (also 
in Table [T]) were fit at the three-loop level to nuclear equilibrium using the 
same method and the Blankenbecler-Sugar propagator for cutoffs of 600 and 
700 MeV respectively. The binding curves for Wl, MOA, MOBl, and M0B2 
are shown in Fig. [61 All four curves are for the total energy up to the three- 
loop level and have been extremized in the both heavy meson fields. One can 
see that while the two- and three-loop level contributions are not large, they 
are not negligible either. Furthermore, nuclear saturation can be reproduced 
at the three-loop level for a range of cutoffs with parameter sets that are 
natural {gs,v ~ 47r). Note that the natural size of the fitted gs implies that 
scalar-exchange diagrams are not redundant. 

Now we want to see how the three-loop energy fits in phenomenologically 
with the power counting. Table [2] shows clearly that the individual terms 
may be large, but strong cancellation occurs. Thus, the proper comparison is 
with the total energy at each order in the loops. Fig. [7] plots the magnitude 
of the energy for both the mean field power counting hierarchy [7.13J along 
side the two- and three-loop energies. Here the three-loop energies utilize the 
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Table 1 

Parameter sets used in this work. Here A is in MeV. 



A 


500 


600 


700 


800 




-20.80 


-28.54 


-35.15 


-40.62 




5.038 


7.658 


9.921 


11.81 


£z-sv{L) 


36.10 


52.73 


68.45 


82.69 


^■i-sviXL) 


- 9.389 


-14.94 


-20.22 


-25.04 


£z-vv{L) 


-16.97 


-27.08 


-38.02 


-49.23 


^■i-vviXL) 


4.438 


7.674 


11.19 


14.80 


£3~Stt{XL) 


0.6105 


1.193 


1.998 


2.997 


£3-Vtt{XL) 


- 1.127 


- 1.544 


- 1.598 


- 1.216 




-5.213 


- 9.503 


-15.05 


-21.84 




1.705 


3.411 


5.617 


8.292 




0.2618 


0.3933 


0.5507 


0.7361 


£3 


- 5.345 


- 8.548 


-12.32 


-16.62 



Table 2 

Individual contributions to the three-loop energy density from ladders (L), crossed 
ladders (XL), and football (F) terms for a range of cutoffs using the Blankenbecler- 
Sugar propagator and the MOA set at nuclear equilibrium [24j . All entries are in 
MeV. 

Blankenbecler-Sugar propagator and cutoffs of 600 and 700 MeV. The total 
two-loop energy is equivalent to third order in the mean field power counting, 
while the three-loop energy is comparable to fourth order. The same hierarchy 
is observed for the Thompson and Erkelenz-Holinde propagators (in fact it is 
more pronounced, as is clear from Fig. E]). Here the mean field parameter sets 
used (Ql and Q2 [7]) include various nonlinearities in the isoscalar meson 
fields; we stress that these nonlinearities were not included in the present 
analysis of the two- and three- loop energies and were shown just to illustrate 
the scale of the mean-field nonlinear terms. 
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Fig. 6. Comparison of the nuclear binding curves for the set Wl (one-loop param- 
eters but with the two- and three-loop contributions included and A = 600), MOA 
(two-loop parameters but with the three-loop contributions included and A = 600), 
MOBl (three-loop parameters with A = 600), and M0B2 (three- loop parameters 
with A = 700). 



It is also of interest to investigate how the individual and total three-loop 
contributions behave with respect to density. Figs.[H]and[9]graph the individual 
ladder integrals and the total three-loop energy as a function of density (p/po 
where p = ps) for the cutoffs 600 and 700 MeV respectively. One can see 
that while the total three-loop energy is roughly linear, the individual ladder 
contributions are nonlinear. 

It is worthy of note that these results disagree with the claim in |44j, that 
the two- and three-loop pion-exchange terms are critical for providing nuclear 
binding and saturation. This follows naturally, since the two-loop pion ex- 
change terms are repulsive, and the (attractive) three-loop terms are smaller 
than the two-loop terms for reasonable values of the cutoff A. The sum of 
the pure pion terms for the MOBl and M0B2 sets is +6.7 and +3.3 MeV 
respectively. 

In summary, we have conducted the loop expansion to third-order for QHD. 
Since covariant many-body loop expansions have been studied for more than 
thirty years [H], it is important to enumerate the new features of this work. 
First, the loop expansion was conducted while incorporating both the chiral 
and heavy meson dynamics. Second, it was shown that, for a reasonable choice 
of cutoff, the energy contributions decrease order by order in loops. Third, 
that natural sets of parameters exist at the three-loop level; note that natural 
parameter sets also exist in the full Brueckner calculation [50j. Fourth, this is 
not a Brueckner calculation, so there was no guarantee the three-loop terms 
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tribution. The abscissa represents the order v in the power counting at the mean 
field level 0. Absolute values are shown. 




Fig. 8. Behavior of ladder integrals as a function of density using the Blankenbe- 
cler-Sugar propagator. Here the MOBl parameter set was used. 



would be small. Fifth, the local contributions appear to be contained in the 
parameterization of QHD. To fully decide whether or not the loop expansion 
is asymptotic for QHD, one needs to consider the ladder and ring sums. 
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cler-Sugar propagator. Here the M0B2 parameter set was used. 
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